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Abstract 
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and another due to the distortion of the Sm-C* . 
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I. INTRODUCTION 



Chiral smectics have two proper characteristics, they are allowed to be ferroelectric and to 
present a helical precession of the optical axes around the layer normal when a tilt of the molecules 
appears in the layers [1]. Like other mesophases tM, they are ferroquadrupolar phases in the sense 
that a large amount of the individual dipoles orient themselves collectively in the bulk and sum 
up in an antiparallel way to give sizeable effects like the flexoelectricity [3]. The liquid crystals 
molecules bear polar links like C=0, N-0, C=N and delocalized electrons (figure \V^), so they 
present a distribution of dipoles all along their skeleton (figure Hh). The Boulder group has shown 
that the molecular dipoles can be approximated without loss of generality by a longitudinal one 
Pl and a transverse one Pt M\ with amplitudes of several debyes. 

Most of the literature in this field has dealt only with the transverse polarization which is at 
the origin of the ferroelectricity [I, sj- 17], and only a few have recognized the importance of the 



longitudinal one [|4|, 



m. 






Figure 1 : dipolar distribution from real molecule to the Boulder model. 



In this paper, we will try to explain the mechanisms at the origin of the formation of the different 
tilted chiral smectics. 



n. THE CHIRAL SMECTIC PHASES 



By order of increasing complexity one encounters the following phases which structure is best 
described by the distorted clock model mainly developed from the data of resonant X-rays scatter- 
ing experiments lll9l - l25ll . 
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A. Sm-A 







Figure 2: Schematic arrangement of molecules in the Sm-A phase. The transverse dipoles vanish while 
the longitudinal ones are in equal numbers up or down. Note that each of the sketched molecules represent 
symbolically one half of all molecules in the bulk, it is not a microscopic view. 



The initial phase which precedes the various tilted phases at higher temperature is the smectic 
A (Sm-A). The molecules are normal to the layers. The transverse dipoles average to zero due 
to a uniform rotation about the long axis. The longitudinal dipoles adopt equiprobable up and 
down orientations (fig. Ill) ensuring that there is no macroscopic polarization but a macroscopic 
uniaxial quadrupole Qij. The uniaxial orientational order parameter (OOP) is expressed as Sij = 
riiTij — where it is the director. When it is written in a frame for which the normal to the 
smectic layers is taken as the z direction it reads [|20 : 



S.. 
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B. Sm-C* 



If the preferred layer thickness decreases with temperature and becomes smaller than the length 
of the molecules, they have to tilt in one direction giving in the simplest case the phase predicted by 
Meyer [ l|, the smectic C* (Sm-C*) where all the molecules are parallel (figure[3]). The transverse 
dipoles give birth to the macroscopic polarization Ps when summed up over at least ten layers. 
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Figure 3: Tilted molecules in a Sm-C* layer, with the longitudinal dipoles in equal number in opposite 
directions while the hindered rotation leaves an average transverse dipole. 



The longitudinal ones have to average to zero but they still sum up in a macroscopic quadrupole 
which main axis is tilted with respect to the layer normal. If one approximates the OOP of the 
Sm-C* as being the same Sij as in the Sm-A tilted at an angle 9 with respect to the layer normal 
in the azimuthal direction $o, one gets : 
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Due to the chirality, the structure precesses around the layer normal z following the law $o = 
qxz with a pitch in the micron range. 

The macroscopic quadrupole ©jj will be to first order proportional to Qij. The more realistic 
case of biaxial is treated in appendix B and keeps the same symmetry as in equ.© with slightly 
involved factorized coefficients. 
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This expression with three basic matrices will be found in all the tilted phases and is funda- 
mental for the continuum theory we have developed. 



C. Sm-C* 




Figure 4: in the anticlinic Sm-C^ phase, the transverse and longitudinal dipoles are compensated when 
pairing the layers, so the phase is ferro-quadrupolar. 



In the anticlinic phase with a period of two layers (figure Hj), both longitudinal and transverse 
dipoles contribute to a macroscopic biaxial quadrupole which has the three C2 directions x, y and 
z as symmetry axes. This phase is misleadingly referred to be antiferroelectric due to the alternate 
orientations of Pt, one should notice that there is also an altemance of Pl that leads to an other 



periodic array in the x direction MliSQ. So it is better to characterize this phase by its quadrupole 
where Pt contributes to Qyy and P^ to ©a.^. and Q^z 

As the Sm-C^ is built by combining <y9 = $o and ip = + tt, the OOP reads : 
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with again a precession $o = (h^ where g2 has the opposite sign to q\ in a given compound. 



D. Sm-C 
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Figure 5: left : side view of the 3 layers unit cell of a Sm-C^j^. at right top view assuming a clockwise 



rotation, the difference in azimuthal angles between layers 1 and 2 or between 2 and 3 is taken as \i ||23 
The in-plane projection of the director of layer 1 makes the angle with x. 



2611 



This phase presents a unit cell commensurate to three layers with unequal changes of the az- 
imuthal angle from layer to layer (A(/9 = or 2(7r — /i)) see e.g figure © with the consequences 
that there is a neat polarization at larger scale (the Sm-C^j^^ is ferrielectric) and a macroscopic pre- 
cession of the structure around the layer normal($o = q'l-^)- AH these informations can be gathered 
when writing the OOP of the phase : 
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taking the definitions of [i and $o given in the figure ([5]), one finds |26[1 that the polarization Ps 
is proportional to / = (1 + 2 cos /u)/3 while the macroscopic quadrupole 6jj is a function of Q, I 
and J = (1 + 2 cos 2/i)/3 and is tilted with respect to the layer normal. 
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E. Sm-C 



Fi2 




Figure 6: left : side view of the 4 layers unit cell of a Sm-C*p^2- at right top view assuming a clockwise 



rotation, the difference in azimuthal angles between layers 1 and 2 or between 3 and 4 is taken as v ||23 , 
The in-plane projection of the bissectrix of layers 1 and 2 makes the angle $0 with x. 



The unit cell is commensurate to four layers with unequal changes of the azimuthal angle from 
layer to layer (Acp = v or n — v) see e.g figure ^ with now no neat polarization at larger scale 
(the Sm-C|^j2 '^ot ferrielectric !) and a macroscopic precession of the structure around the layer 
normal ($0 = ^^2^^ almost all the studied compounds, qi and q2 have the same sign while q2 and 
q'l have the opposite). All these informations are gathered in the OOP of the phase : 
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with the definitions of v and $0 given in the figure Q, one finds [|26h that the macroscopic 
quadrupole Ojj which has the layer normal as one of its eigenaxes is a function of 9 and 

J = — cosv. 
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F. Sm-C 



d6 



A last commensurate phase with six layers has been predicted by H&T [1231 . 12711 and recently 
evidenced by Shun Wan et al [12 811 . We will not develop on it but it has a symmetry close to that of 
Sm-Cpj2 similar properties. 



G. Sm-C! 



Last but not least, this phase shows a periodic precession with a short period which is not 
commensurate to the layer thickness. Its macroscopic OOP is simply uniaxial : 
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H&T have shown that the Sm-C* is fundamental for the obtention of the commensurate subphases. 
In short the incommensurate period varies continuously with the temperature taking values com- 
prised in a subset of the interval of 2 to 8 layers. They have shown [|23L 12711 that when this period 
gets close to an integer number of layers, the system prefers to lock-in at this integer value at the 
expense of the twist energy compensated by some other gains. They have proposed that this gain 
scales at due to an anisotropy of in-plane elastic energy. What we propose here is that the lock- 
in allows the onset of macroscopic quadrupole and sometimes dipole with and P contributions 
that explain the development of the full set of subphases. 



III. THE SEQUENCE OF TILTED PHASE 



A. The macroscopic and microscopic orientational order parameters (OOP & oop) 

We have just seen that all the phases described by the distorted clock model are characterized by 
their macroscopic orientational order parameter (OOP) Qij which general form valid in all phases 
has been given in equation dH). It is defined on a scale of at least ten layers like the quadrupole 
Qij and the polarization i^. We have also recalled the fundamental statement of H&T that the 
helicity of the Sm-C* governs the appearance of other phases. So one has to develop the theory of 
the Sm-C* phase and its transition from the Sm-A. For that we consider that each layer in a tilted 



8 



smectic phase is such that the director makes an angle 6 with z while its in-plane projection makes 
the angle (p with a; ; so we express the result of the rotation of Sij in the xyz frame as the tensor 
Sij which is the microscopic orientational order parameter (oop) of the layer : 
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As already reported in [|29, [s^ the oop splits into three traceless invariants, namely a bulk 
3D-uniaxial tensor which depends only on 9, and two in-plane tensors respectively 2D-biaxial and 
2D-uniaxial depending also on tp. This local tensor will be used to compute the bulk OOP of each 
tilted phase by including the z dependence of the azimuth angle cp. 



B. Landau - de Gennes free energy 

The free energy density describing the phase transition from the Sm-A to tilted phases can be 
written as a power series of the local oop Sij : 
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Following the Smith and Rivlin theorem 13 1|] we express the tensorial coefficients aiju, Vlijkimn 
and bijkimnop as products of the elementary tensors like the Kronecker 5ij, the Sm-A OOP Sij, the 
vacuum tensor Vij and as we deal with chiral compounds the fully antisymmetric Levi-Civita odd 
tensor eijk- 



After some tedious calculations |2^ 3^ 
the invariants introduced in (|7]) : 



one gets rather simple results which are functions of 
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using the approximation sin 6* oc 6* in (|7]), the free energy density can then be developed in a 
power series in 9, without any dependence on the azimuth This is not surprising as 9 and can 
be considered as the modulus and the phase of the Sm-C complex order parameter 9 exp{i(p) [|2|] : 



F, = ^a^9^ + h)9^ + ... (12) 

where as = 03 (T — T^) governs the Sm-A to Sm-C phase transition. This 9'^ term comes from 
the 2D-uniaxial invariant while the 9^ one is the sum of ai, 02, i^i, ^^3 and 63 contributions. 



In the mean-field approximation, the angle 9 behaves like ^ a/Tc — T in the tilted smectic 



phases. This looks like previous theories Da, 



l3,[l5l 



which make use only of the 2D-uniaxial 



invariant missing somewhat the other terms which are at the origin of the Hamaneh-Taylor (H&T) 



theory [|23 



From now on we have considered only the case of constant azimuthal angle tp like in achiral 
smectics. As we are dealing with chiral compounds, we can introduce the gradients of (/9 in the z 
direction to take into account the helicity of chiral smectics. 

C. helicity of the first tilted phase 

When the tilt appears at the transition, the first phase that condenses can be the Sm-C* or the 
Sm-C* , they are distinguished by the value of the pitch of the helix. We have to add to the Fi term 
the energy contributions of coming from the gradients of the oop : 
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The only gradients which do not perturb the symmetry of the layers are the twisting ones which 
induce a rotation of the 2D-uniaxial and biaxial invariants around the layer normal. So we assume 
that the first term linear in measures the twisting power of the phase while the second one, 
quadratic, reflects the cost in twist elastic energy. One finally gets : 



^ijklm ^ij^ k^lm -^1 {.i^^xx ^yy^dzSxy ^xy^zi^^xx '^J/J/)) ~^ ^2i^^xzdzSyz Sy^d^Sxz^ (14) 
kijklmrNiSjk^lSran = h {{d^Sccx " d^Syy)^ + 4 (d^SxyY) + ^2 {{d^S^^Y + {dzSy^f) (15) 

when reduced to functions of 9 and ip it remains : 

F^ = F^-\d-,^ + h{d,^f (16) 

with the trivial solution ip = (pQ + qz describing an helicoidal rotation at the wave vector q = \/k. 
The introduction of helicity renormalizes slightly the quadratic term in the free energy leading 



to a small displacement of the transition temperature Tc to Td [11711 . 

Two cases must be distinguished now, depending on the magnitude of the helical pitch p = 

2n/q. 

1. small pitch : the Sm-A to Sm-C*^ phase transition 

When the helical pitch takes values typically from 2 to 8 layers, the Sm-C* phase appears at the 
transition |32]. Although it has been reported in a few compounds a Sm-C* phase with 15 to 50 
layers 11231 13211. some caution must exercised as it is on the basis of D.S.C. data taken at 3° C mn^^, 
with a small bump appearing above a large peak, that the phase has been reported [33]. To be sure 
of the phase identification this bump should have been followed at lower speeds as was done in 



MHPOBC 13^ 

The Sm-C* phase is fundamental in order to get the full sequence of commensurate subphases. 
It governs the number of layers of the unit cell by means of the azimuth increment a as developed 
in H&T theory. One of the experimental challenges risen by this theory is the measurement of 
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a = 271 d/p by means of the pitch p of the Sm-C* phase. Let us point out that q = 2'7i/p is the 
ratio of two polynomials in 6 without any critical dependence : 



we expect it to vary smoothly from \2lk2 at the transition to finite values later on. In all the 



known experiments to date lll9 . 



21 



22|, 



32|, 



35 



36(1 there is a general trend of increasing 



a 



qd ~ qdo(l — 9"^ /2) when cooling down at the exception of one compound |13 



3311 where 



the Sm-C* denomination subject to caution. Let us remark that the thermal variation of qd is the 
product of two terms and cannot be predicted for sure. Close to the transition, the decrease of d 
when cooling down may be dominant but it seems that later on q will increase and at the end will 
determine the response. 



2. large pitch : the direct Sm-A to Sm-C* phase transition 



When the preferred pitch is typically larger than 0.3 fim, the tilted phase can be considered 
locally as a Sm-C* which precesses slowly around the layer normal. One has then to take into 
account the macroscopic polarization [1] which rotates too. In order to be coherent, P is defined 
over a few layers (~ 10) , while Sjk is relative to one layer and the layer polarization is not a 
macroscopic quantity pTI] . One has to introduce the macroscopic OOP Qjk =< sjk > which is an 
average over the same area. This ensures that when the pitch is smaller than 10 layers, i.e. in the 
Sm-C* phase, both P and the extra parts in < sjk > vanish. As already stated earlier Is, 29, 3o|] . 



the polarization can be formally introduced within linear couplings with the 2D-uniaxial invariant 
of the OOP, one describing the ferroelectricity [1] and the other the flexoelectricity Is, 38] : 
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The total polarization is obtained by the minimization with respect to P of the following en- 
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ergy : 



One gets : 
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As expected by symmetry the polarization is an in-plane vector which maximizes AFp when 
it is normal to the projection of the director (Qxz, Qyz) or equivalently parallel to the gradient 
{dzQxz, dzQyz)- Both contributions to P are coUinear, they follow the helical precession of the 
director and they change sign with the chirality. The polarization which is measured usually in 
unwound samples is the ferroelectric one, as the other contribution disappears 117(1 . 

Here again the introduction of the macroscopic ferroelectric polarization renormalizes slightly 
the quadratic term in the free energy leading to another small displacement of the transition temper- 
ature Tci to OOn . Conversely, the flexoelelectric polarization changes the twist elastic constant 
and the helical pitch. 



IV. THE COMMENSURATE SUBPHASES 



The Sm-C* and Sm-C* phases are not the only ones encountered in these compounds, when 
further cooling down a sequence of commensurate phases with unit cells of 1 to 6 layers have been 

nn o 

reported 151,1280 which are best described by the distorted clock model [|19|] . The fundamental idea 
in H&T theory is that the basic tilted phase obtained below the Sm-A phase is the Sm-C* one, with 
a short pitch varying from about 2 to 8 layers. When the value of the pitch is close to an integer 
number of layers, there can be a lock-in of the structure at this integer number at the expense of 
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the twist energy, provided that there is a gain in electrostatic or elastic energy [|23L 

Let us enforce the fact that in these subphases, the tilt angle 6 and the a parameter are functions 
of the temperature only given by the resolution of equations (fT2l) and (fT6l) respectively. 
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Figure 7: Sketch of azimuthal angles in the unit cells of Sm-C^-^ (top) and Sm-C^-2 phases (bottom). The 
real distribution in the center is a combination of XY (left) and Ising models (right). 

A. Hamaneh-Taylor theory (H&T) 



The commensurate subphases of the distorted clock model are a compromise between thepure 



m 



XY clock model 111911 that requires a regular increase of the azimuth and the old Ising model O, 
that obliges the tilt to lie in one direction of the plane (see e.g. figure ([7])). The azimuthal angles in a 
subphase deviate from the regular XY behavior at the expense of the twist elastic energy while the 
Ising model favours in-plane anisotropy. H&T have introduced the angles A(y9 and a to describe 
the distance from the XY model and the parameter J = (cos 2^9) < 1 to describe the biaxial 
Ising-like tendency 11231 12711. We have further defined the parameter — 1 < / = (cos v^) < 1 that 
measures the in-plane uniaxial character 

I and J have already been introduced in the definition of the OOP Qij =< Sij > in equation 
(HI), which is a function of the tilt angle 9 and of $0 defined as the angle between the origin of 
azimuthal angles in the unit cell and the x axis ll26ll . The resulting order parameter Qy is unique ; 
it is only its expression in a given frame which depends on $0. 



B. cost of lock-in : the short range term 



From equation (1161) one knows the amount of twist energy per unit volume lost when the wave 
vector is slightly different from its preferred value q = 27r/p. It amounts to AFt = k{d(f/dz — 
qY /2. Introducing the azimuthal angle increment Aip and the layer thickness d, one gets AFt = 
k{Aip — qdY 126?, AFt has a zero minimum value when A(/9 = qd = a. 
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This expression averaged over the unit cell of any distorted clock model structure can be related 



to the short range term Fsr introduced empirically by H&T [|23ll . Fsr = Fq < cos(A(/? — a) > 
is a positive energy which takes its maximum value Fq in the Sm-C* phase when A^p = a. On 
taking Fq = k/dP and cos(A(y9 - a) ~ 1 - {Aip - of jl one gets Fsr = Fq - AFt. 

We then know how much it costs to lock- in the average increment < A^p >= 27r/n at a value 
close to a. In the initial clock model this increment was supposed to be constant in the locked-in 
phases but they would only be commensurate Sm-C* phases without any change in energy. In that 
case the unit cells are uniaxial without in-plane anisotropy so that the extra parameters I and J are 
identically null. 

C. gain from quadrupolar and dipolar ordering : the P and ,P terms 

Uniaxial nematic and Sm-A phases are well known to be ferro-quadrupolar [3] i.e. the 
molecules are arranged at rest so that their microscopic electric dipoles and quadrupoles sum 
up cooperatively to give a macroscopic quaqrupole density proportional to the uniaxial OOP of 
the phase : 

e° = s,, (22) 

This can be demonstrated the following way : the quadrupolar density has a quadratic self- 
energy 0^^/2x0 and is linearly coupled to the OOP by the term —<da'S>ijSij/2xe^ this leads to 
the result given in equation (l22l) and to the expression of the energy gained by the creation of the 
quadrupolar density : 

^Fi = (23) 
In the Sm-C* phase, this energy decreases with 6 as : 



AFg = ( 1 - ^ sin^ e ] (24) 



3xe V 2 

the lock-in to commensurate subphases allows to recover a part of this loss. In the uniaxial ap- 
proximation we are using, the quadrupolar energy in the Sm-C* phase is the same as in the Sm-A, 

AF^ = AFi so : 

AF^ = (j^l-^sin^e^ + ^sin^^ + Ssin'^cos'e 1 (25) 
a straightforward generalization to any value of / and J reads : 
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We are eventually left with the gain of lock- in which reads : 
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When completed by the similar term due to the presence of the macroscopic polarization Ps if 
/ 7^ as developed by Dhaouadi et al. [I26], one gets : 
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the full energy gain reads with this n ew P term : 



(28) 
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D. balance between short range loss and long range gain 

When a phase described by the distorted clock model appears, it is characterized by non zero 
values of / and J that minimize the following energy at a negative value : 

'1 



{Aif - ay - r]J' - ^^I' 



(30) 



The phase diagrams in the plane (0 ^7r,0^77^1) have been computed by H&T 



[123 



2711 with the term and Dhaouadi [26] (/^) together with their behavior under an applied 



electric field {-PE) in the last case. 

Let us remark that up to now we have discussed the commensurate subphases in the unwound 
geometry although we know that they are all precessing around the layer normal. We propose to 
treat this problem now and compare our results with the well known pitch inversion at the Sm-C^^j^ 
to Sm-Cpj2 phase transition. 
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Figure 8: Sketch of the short range term expressing the cost of lock-in from Sm-C* (A(/? = e) to Sm-C* 
phase (A(^ = 0). A small positive increment 6{Aip) would reduce the cost. 

E. helicity of the subphases 

We have considered when computing the phase diagrams that the unit cells of the commen- 
surate subphases were frozen and did not rotate while the subjacent Sm-C* phase experiences a 
strong spontaneous twist. We first examine graphically the simple cases of the Sm-C* and Sm-C^ 
phases where the azimuth increase between layers is equal to Aip = or Aip = ir while the Sm-C* 
phase that would take place otherwise has a small value a = e in the first case (figure [8]) or is close 
to vr in the second (a = n — e figure ID). The figures show that the cost in twist energy of the lock- in 
will be reduced if the Sm-C* rotates in the same direction as the Sm-C* while the Sm-C^ has to 
take the opposite sense. This illustrates the general trend observed in the experiments [S] that the 
sense of the helix is opposite in the two phases for a given compound. We have thus shown that 
the would-be Sm-C* phase exercises a kind of torque on the azimuthal angle $o with a non trivial 
sign. It has to be completed by a spontaneous twist we will develop now. 

We first take into account the preceding remarks by replacing AFt = Fq (Ay? — a)^/2 by 
AFt = Fq {Aip + d (9^$o — af 12. We then state that the macroscopic angle $o is subject to the 
same laws than the microscopic : 

Ft = Fo {A^ + d 9,$o - a) 72 - A + ]^K{d,^^f (31) 
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Figure 9: Sketch of the short range term expressing the cost of lock-in from Sm-C* (A99 = tt — e) to Sm-C^ 
phase {Aip = tt). A smaller Atp would reduce the cost. 



replacing F^hy k/(P, one has to minimize : 



yielding : 



-(k + K) (9,$c 



(A<^ - a) - A d, $0 
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thus the wave vector of the macroscopic helicity is given as usual by the ratio of the twist sources 
to the rigidity. We already know that the pitches are much larger than in the Sm-C* phase so one 
may assume safely K k. The source in the numerator is made of two terms. The first one is 
the intrinsic twisting power A which we take as having the same sign as the microscopic one A, 
this remaining to be checked in real world. The second one corresponds to the discussion we just 
developed with the help of figures ([8]) and 

Let us briefly determine the sign of dz^o in the different commensurate subphases : 

1. in Sm-C* the second term dominates and has the sign of a. 

2. in Sm-C^ the second term dominates with the opposite sign. 



3. in Sm-C|;.j]^ if the angle yU is much larger than 27r/3 as reported in literature [|26l] . the opposite 
sign comes again. 
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4. the Sm-C^j2 more involved as it seems that the first source dominates and in our hypoth- 
esis on the sign of A, the helix has the same sign as a. 

So we have found that usual sequence with a sign change in the middle of the range is well 
explained in our theory. 

V. CONCLUSION 

We have demonstrated that the sequence of distorted commensurate phases observed in tilted 
chiral smectics is explained by the gain in electrostatic energy due to the lock-in of the unit cell 
to a number of layers which is the integer closest to the ratio pitch over thickness of the subjacent 
Sm-C* phase. We also explain the sign change of the helicity in the middle of the sequence by a 
balance between two twist sources. 

Appendix A: shifts 

1. due to the helicity 

remembering equation (fT6l) and introducing the preferred wave vector q = \/k, one gets F2 = 
Fi — kq^ /2 = F\ — }? I2k. Fi is an even polynomial in Q as well as the correction term, the 
leading term reads now after equation (O a^iT - Tc)Q'^l'l - \\l'lk2 = a^{T - T*)e'^/2 with a 
straightforward definition of T* > T^. 

2. due to the polarization 

Here again the correction term due to the macroscopic ferroelectric polarization can be ex- 
pressed as AFp = —xC^O'^ leading to another shift in the transition temperature AT* = xC^/ois. 
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Appendix B: biaxiality vs uniaxiality 



All the unwound phases in the distorted clock model are biaxial with the in-plane C2 axis as 
one eigenaxis. Their OOP Sy and the rotation matrices should read : 



V 



-a 
0-6 
a + b 



\ 



( 



cos 6 sin^^ 

1 
- sin 6 cos 6 



\ 



^ cos (p sin (y9 ^ 



(Bl) 



Pij = 10 = —simp cos (fi 



when one executes a rotation of angle 9 around this axis and another of angle ip around 3, one 
gets : 



Q 



- (a cos 26 + b cos^ 9) 



1 



■1/3 ^ 
-1/3 
+2/3 J 



^ cos 2ip sin 2(p ^ 



+ ^{-acos2e + b{l + sm^e)) 



( 



— (2a + b) sin Q cos d 



sin 2^9 — cos 2(y9 


\ 



(B2) 



cos ip 
sin 
y cos sin I 

this differs only slightly in the coefficients from the uniaxial form we have used. We took 
a = b = 1/3, we could have been closer to reality with a = 1/3, b = 1/3 + 5 because of Lf 
fluctuations evidenced by conoscopy under field |15|]. 



Acknowledgments 



We wish to acknowledge the support of CRPP's pot commun and the help of the Tunis group 
for this work. 



20 



[1] R. B. .Meyer, L. Liebert, L. Strzelecki, and P. Keller, J. Phys. (Paris) 36, L69 (1975). 
[2] P. G. de Gennes and J. Prost, The Physics of Liquid Crystals (Clarendon Press, Oxford, 1993), 2nd ed. 
[3] J. Prost and J. P Marcerou, J. de Phys. (Paris) 38, 315 (1977). 
[4] D. R. Link, J. E. MacLennan, and N. A. Clark, Phys. Rev. Lett. 77, 2237 (1996). 
[5] A. D. L. Chandani, E. Gorecka, Y. Ouchi, H. Takezoe, and A. Fukuda, Jpn. J. Appl. Phys. 28, LI 265 
(1989). 

[6] J. Lee, A. D. L. Chandani, K. Itoh, Y. Ouchi, H. Takezoe, and A. Fukuda, Jpn. J. Appl. Phys. 29(6), 
1122(1990). 

[7] A. Fukuda, Y. Takanishi, T. Isozaki, K. Ishikawa, and H. Takezoe, J. Mater Chem. 4, 997 (1994). 
[8] N. M. Shtykov, A. D. L. Chandani, A. V. Emelyanenko, A. Fukuda, and J. K. Vij, Phys. Rev. E 71, 
021711 (2005). 

[9] T. Carlsson, B. Zeks, A. Levstik, C. Filipic, 1. Levstik, and R. Blinc, Phys. Rev. A 36, 1484 (1987). 
[10] H. Orihara and Y. Ishibashi, Jpn. J. Appl. Phys. 29, LI 15 (1990). 
[11] H. Sun, H. Orihara, and Y. Ishibashi, Jpn. J. Appl. Phys. 60, 1991 (1991). 
[12] A. Roy and N. Madhusudana, Europhys. Lett. 36, 221 (1996). 
[13] A. Roy and N. Madhusudana, Eur. Phys. J. E. 1, 319 (2000). 
[14] N. Vaupotic and M. Cepic, Phys. Rev. E 71, 041701 (2005). 
[15] V. Lorman, Liq. Cryst. 20, 267 (1996). 

[16] H. F. Gleeson, L. Baylis, W. K. Robinson, J. T. Mills, J. W. Goodby, A. Seed, M. Hird, P. Styring, 

C. Rosenblatt, and S. Zhang, Liq. Cryst. 26, 1415 (1999). 
[17] G. Durand and P. Martinot-Lagarde, Ferroelectrics 24, 89 (1980). 

[18] T Matsumoto, A. Fukuda, M. Johno, Y. Motoyama, T. Isozaki, T. Yui, S. S. Seomun, and M. Ya- 

mashita, J. Mater Chem. 9, 2051 (1999). 
[19] P Mach, R. Pindak, A. M. Levelut, P. Barois, H. T. Nguyen, C. C. Huang, and L. FurenUd, Phys. Rev. 

Lett. 81, 1015 (1998). 

[20] A. Cady, J. A. Pitney, R. Pindak, L. S. Matkin, S. J. Watson, H. F. Gleeson, P. Cluzeau, P. Barois, 

A. M. Levelut, W. Caliebe, et al, Phys. Rev E 64, 050702(R) (2001). 
[21] L. S. Hirst, S. Watson, H. F. Gleeson, P. Cluzeau, P Barois, R. Pindak, J. Pitney, A. Cady, P M. 

Johnson, C. C. Huang, et al., Phys. Rev. E 65, 041705 (2002). 



21 



[22] D. A. Olson, X. F. Han, A. Cady, and C. C. Huang, Phys. Rev. E 66, 021702 (2002). 
[23] M. B. Hamaneh and P. L. Taylor, Phys. Rev. Lett. 93, 167801 (2004). 

[24] Z. Q. Liu, B. K. McCoy, S. T. Wang, R. Pindak, W. Caliebe, P. Barois, P. Femandes, H. T. Nguyen, 

C. S. Hsu, S. Wang, et al., Phys. Rev. Lett. 99, 077802 (2007). 
[25] N. W. Roberts, S. Jaradat, L. S. Hirst, M. S. Thurlow, Y. Wang, S. T. Wang, Z. Q. Liu, C. C. Huang, 

J. Bai, R. Pindak, et al., Europhys. Lett. 72, 976 (2005). 
[26] H. Dhaouadi, N. Bitri, S. Essid, T. Soltani, A. Gharbi, and J. P. Marcerou, Phys. Rev. E 80, 031712 

(2009) . 

[27] M. B. Hamaneh and P. L. Taylor, Phys. Rev. E 72, 021706 (2005). 

[28] S. Wang, L. Pan, R. Pindak, Z. Q. Liu, H. T. Nguyen, and C. C. Huang, Phys. Rev. Lett. 104, 027801 

(2010) . 

[29] J. R. Lalanne, J. Buchert, C. Destrade, H. T. Nguyen, and J. P. Marcerou, Phys. Rev. Lett. 62, 3046 
(1989). 

[30] L. Dupont, M. Glogarova, J. P. Marcerou, H. T. Nguyen, C. Destrade, and L. Lejcek, J. Phys. 11 

(France) 1, 831 (1991). 
[31] G. F. Smith and R. S. RivUn, Ns- DA- 3487/3, Brown University (1955). 

[32] A. Cady, D. A. Olson, X. F Han, H. T. Nguyen, and C. C. Huang, Phys. Rev E 65, 030701 (R) (2002). 

[33] V. Faye, J. C. Rouillon, C. Destrade, and H. T Nguyen, Liq. Cryst. 19, 47 (1995). 

[34] J. Philip, J. R. Lalanne, J. P Marcerou, and G. Sigaud, Phys. Rev. E 52, 1846 (1999). 

[35] A. Cady, X. F. Han, D. A. Olson, H. Orihara, and C. C. Huang, Phys. Rev. Lett. 91, 125502 (2003). 

[36] V. Laux, N. Isaert, G. Joly, and H. T. Nguyen, Liq. Cryst. 26, 361 (1999). 

[37] J. P Bedel, J. C. Rouillon, J. P. Marcerou, H. T. Nguyen, and M. F. Achard, Phys. Rev E 69, 061702 

(2004). 

[38] R. B. Meyer, Phys. Rev. Lett. 22, 918 (1969). 

[39] F. Beaubois, J. P. Marcerou, H. T. Nguyen, and J. C. Rouillon, Eur. Phys. J. 3, 273 (2000). 



22 



